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Abstract 

We reanalyse the time-variable lightcurves of the transiting planetary system PTFO 8-8695, in which a 
planet of 3 to 4 Jupiter mass orbits around a rapidly rotating pre-main-sequence star. Both the planetary 
orbital period P or b of 0.448 days and the stellar spin period P s less than 0.671 days are unusually short, 
which makes PTFO 8-8695 an ideal system to check the model of gravity darkening and nodal precession. 
While the previous analysis of PTFO 8-8695 assumed that the stellar spin and planetary orbital periods 
are the same, we extend the analysis by discarding the spin-orbit synchronous condition, and find three 
different classes of solutions roughly corresponding to the nodal precession periods of 199±16, 475±21, and 
827 ±53 days that reproduce the transit lightcurves observed in 2009 and 2010. We compare the predicted 
lightcurves of the three solutions against the photometry data of a few percent accuracy obtained at 
Koyama Astronomical Observatory in 2014 and 2015, and find that the solution with the precession period 
of 199 ± 16 days is preferred even though preliminary. Future prospect and implications to other transiting 
systems are briefly discussed. 
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1. Introduction 

PTFO 8-8695 is a pre-main-sequence star (weak-line T- 
Tauri star) located in the Orion OBla star forming re¬ 
gion, first discovered by Briceno et al. (2005). Subsequent 
photometric transit and spectroscopic radial velocity mea¬ 
surements revealed that PTFO 8-8695 harbors one close- 
in planet (PTFO 8-8695 b) with planetary mass M p < 
5.5Mj (van Eyken et ah, 2012). This is the first transit¬ 
ing exoplanet candidate around a pre-main-sequence star. 

van Eyken et al. (2012) discovered that PTFO 8-8695 in 
2009 and 2010 exhibited very different transit lightcurves. 
Barnes et al. (2013) (hereafter, B13) found that the pecu¬ 
liar time-variation can be explained by the gravity dark¬ 
ening of the central star and the spin-orbit nodal pre¬ 
cession of the system. Indeed they are successful in fit¬ 
ting simultaneously the lightcurves in 2009 and 2010, and 
estimated the system parameters including stellar and 
planetary radii (i? s ~1.0 Rq, R p ~1.7 i?j), planetary mass 
(M p = 3.0 or 3.6 Mj), and spin-orbit angle (</>~70°; the 
angle between stellar spin and planetary orbital axes). 

In doing so, they assumed that the stellar spin and 
planetary orbital motion are synchronised for simplicity. 
Given the large spin-orbit misalignment estimated for the 
PTFO 8-8695 system, however, the synchronous assump¬ 
tion may not be justified from the dynamical point-of- 
view. Therefore we reanalyse the PTFO 8-8695 transit 
photometry without adopting the spin-orbit synchronous 
assumption. As expected, the permitted parameter space 


is increased, and we discuss a feasibility to break the de¬ 
generacy among different solutions with future observa¬ 
tions. 

The rest of the present paper is organized as follows. 
Section 2 briefly summarizes the previous work on the 
PTFO 8-8695 system. We adopt a model of gravity dark¬ 
ening and spin-orbit nodal precession of B13, which is 
described in Section 3. Our data analysis and parameter 
fitting procedures are presented in §4.1 and 4.2, respec¬ 
tively, and the resulting constraint on the system parame¬ 
ters from 2009 and 2010 data is shown in §4.3. We also dis¬ 
cuss how the long-term photometric monitoring of the sys¬ 
tem is useful in further constraining the parameter space. 
A preliminary analysis with the observation performed in 
2014-2015 at Koyama Astronomical Observatory (KAO) 
is presented in §4.4. We discuss the long-term dynam¬ 
ical stability of the system in Section 5. Final section 
is devoted to the summary of the paper and the further 
discussion. 

2. Exoplanetary system PTFO 8-8695 

The transits of PTFO 8-8695 b were observed at 
Palomar Observatory for 2009 December 1 - 2010 January 
15 (hereafter, referred to as 2009 observation), and 2010 
December 8-17 (2010 observation) as part of the Palomar 
Transient Factory (PTF) that ran for 2009 - 2012. The 
2009 and 2010 observations comprise 11 and 6 reliable 
transit lightcurves, respectively, with the transit period of 
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0.448413 days (van Eyken et al., 2012). As already men¬ 
tioned in Introduction, the lightcurves exhibit a signifi¬ 
cant time variation, and B13 solved this puzzle by taking 
into account both the rapid rotation of the star and short 
transit period of the planet. 

Spectroscopic measurement of PTFO 8-8695 implies 
that the host star is a rapid rotator with the spin period of 
P s < 0.671 days (van Eyken et al., 2012). This is consistent 
with the observational evidence that T-Tauri stars (espe¬ 
cially, weak-line T-Tauri stars like PTFO 8-8695) tend to 
rotate more rapidly than main sequence stars (Edwards 
et al., 1993). Because of the rapid rotation, the equator 
of PTFO 8-8695 would be significantly expanded relative 
to its polar radius. 

Because of the stellar surface distortion, PTFO 8-8695 
should exhibit significant gravity darkening (von Zeipel 
1924, Lucy 1967) and the nodal precession of stellar spin 
and planetary orbital axes. B13 took account of both 
effects properly, and are successful in reproducing both 
2009 and 2010 lightcurves either for two different stellar 
masses, M s = 0.34 and 0.44 M 0 (Briceno et al., 2005). 

The nodal precession period depends sensitively on both 
the stellar spin and planetary orbital periods, P s and P 0 rb- 
While the latter is well-determined from the transit period 
(0.448 days), only the upper limit is obtained observation- 
ally for the former, P s < 0.671 days. Since van Eyken et al. 
(2012) concluded that the stellar spin and planetary or¬ 
bital motion are most likely synchronized (P s = P or b), B13 
adopted that assumption for simplicity throughout their 
analysis. 

Theoretically speaking, however, the validity of this as¬ 
sumption is not clear. One of the most probable mech¬ 
anisms for the spin-orbit synchronization is the tidal ef¬ 
fect between the star and planet. In this case, the spin- 
orbit synchronization and the spin-orbit alignment oc¬ 
cur simultaneously. The equilibrium tidal model, for in¬ 
stance, predicts that the effective timescales for the two 
processes are almost the same (Lai 2012, Rogers & Lin 
2013). Nevertheless the result of B13 indicates the signif¬ 
icant spin-orbit misalignment (4>~70°), and indeed they 
admitted that the truly synchronous rotation might be 
difficult to be achieved under such a misaligned state. 
This is why we attempt in the present paper to reanal¬ 
yse the PTFO 8-8695 system, to see how the estimated 
parameters are sensitive to the synchronous condition. 

3. Basic equations 

This section summarizes the basic equations that we use 
in the later analysis, largely following B13. Our model 
employs 16 parameters in total, and for clarity we list 
them in Table 1. 

3.1. Configuration of the system 

Figure 1 defines the parameters specifying the geometric 
configuration of the system. We define the Cartesian coor¬ 
dinates in such a way that the star is located at the origin 
and the x — y plane coincides with the sky plane with the 
positive 2 -direction pointing towards the observer. The 


stellar spin vector projected onto the sky plane is defined 
as the positive y-direction, and the ^-direction is set so as 
to form a right-handed triad. 

The stellar inclination and planetary orbital inclination, 
i s and i OI b, respectively, are measured from the line-of- 
sight pointing to the observer (positive 2 -axis). Longitude 
of the ascending node, f2, is measured from the positive 
z-axis to the direction towards the ascending node coun¬ 
terclockwise. 



Fig. 1. Schematic illustration for spin-orbit angle cp as 
a function of stellar inclination i s , planetary orbital in¬ 
clination i or b, and longitude of the ascending node Q. 

In what follows, S and L denote the stellar spin 
and planetary orbital angular momentum vectors, respec¬ 
tively. We define the total angular momentum vector as 
J = S + L. Their unit vectors and norms are denoted as 
s, l , and j, and S, L, and J, respectively. 

In terms of the angles defined in Figure 1, s and l are 


written as 



( ° \ 

/ 

sin* or b sinfl 

s = I sinig I 

- *= 

— sini or b cosfl 

\ cos z s / 

\ 

COS 7-orb 


Also the spin-orbit angle f> is written as 
f = cos 

= cos _1 (— sun s sini or b cos + cosi s cosj or b)- ( 2 ) 

3.2. Flux profile of the central star 

We compute the stellar flux by numerically integrating 
the stellar intensity profile B\(T)xI(p) throughout the 
visible stellar disk, where B\(T ) is the Planck function 
at the wavelength of A, and I(p) is the limb-darkening 
function. The following analysis is performed at the 
monochromatic wavelength of A = 0.658 pm (see section 
4). We adopt quadratic limb-darkening law in our analy¬ 
sis: 

I(p) = 1 - iti(l - p) - u 2 { 1 - p) 2 , 


( 3 ) 
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Table 1. Parameters for the analysis. 


parameter 

symbol 

fixed value 

note 

Mean stellar density 

Ps 

- 

a 

Stellar effective temperature at the pole 

T P oi 

3470 K 

b 

Stellar rotation period 

P S 

- 


Time of inferior conjunction 

tc 

2455543.9402 HJD 

c, d 

Stellar inclination at t c 

is 

- 

e 

Stellar moment of inertia coefficient 

c 

0.059 

f 

limb-darkening parameter 

Cl = Ul + U2 

0.735 

f 

limb-darkening parameter 

C2=U\— U2 

0.0 


Gravity-darkening parameter 

P 

0.25 

f 

Planet-to-star mass ratio 

M p /M s 

- 


Planet-to-star radius ratio 

Rp/R s 

- 

a 

Planetary orbital period 

P 3rb 

0.448413 days 

d 

Orbital eccentricity 

ecosu; 

0 

g 

Orbital eccentricity 

esino; 

0 

g 

Planetary orbital inclination at t c 

^orb 

- 

e 

Longitude of the ascending node at t c 

n 

- 

e 


a specified with stellar equatorial radius R s ,e q5 M s /^7rilf jeq or ^p/^s,eq- 

b Briceno et al. (2005) reported 3470 K as the stellar effective temperature, but we use that value 


as polar temperature assuming that their difference in the analysis is negligible. 

c an epoch when uj + / = 7r/2. 

d van Eyken et al. (2012). 

e Barnes et al. (2013) adopt different notations for these three angular parameters (stellar obliquity 
ip, planetary orbital inclination i, projected spin-orbit angle A) from ours, and they are related as 
Ip = is — 7T /2, i — 7T — 2 or b i an d A = 7T — f2. 

f Barnes et al. (2013). 

g We assume a circular orbit following Barnes et al. (2013). This assumption is supported by an 
equilibrium tidal theory, which predicts that a close-in planet acquires the circular orbit on a much 
shorter time scale than those for spin-orbit synchronization or alignment. 


where c\ = u\-\-U 2 and C 2 = u± — U 2 are quadratic limb- 
darkening coefficients, and g denotes the cosine of the 
angle between line of sight and the stellar surface normal. 
The effective temperature T at each point of the stellar 
disk is estimated with von Zeipel’s theorem (von Zeipel, 
1924), which says that the temperature distribution on the 
stellar surface induced by gravity darkening is specified by 
the stellar surface gravity g : 

JL = (jl 

Tpoi \ 3poi 

where T po i and g po \ are the effective temperature and sur¬ 
face gravity at the pole, and /3 is the gravity-darkening 
parameter. See Barnes (2009), Barnes et al. (2011), and 
Masuda (2015) for more details of our model. 

3.3. Spin-orbit nodal precession 



We solve the orbit-averaged equations of motion for the 
two rigidly-rotating bodies (Boue & Laskar, 2009): 


S = ~lX8, (5) 

i = — ysxl, (6) 


where 

3 GM P 

P ~ 2a 3 (l-e 2 ) 3 /2 


3 2 M S M P 1 

— TL --- 

2 M s + M p (1 — e 2 ) 3 / 2 

R s \ 2 1 

a 


■■ ncoscj)-J 2 


(1 — e 2 ) 2 


J 2 RI COS (j) 


xL 


: —flpL. 


(7) 


Here a, n = A JG(M S + M p )/a 3 = 27 r/P or b, and J 2 denote 
orbital semi major axis, orbital mean motion, and the 
second-order gravitational coefficient of the central star, 
respectively. In this formulation the planet is considered 
as a point mass (i.e., planetary spin is neglected) and 
higher-order gravitational moments are neglected. 

Equations (5) and (6) are rewritten in terms of J = 
S + L = Ss + LI: 


s = -^r(J-Ss)xs 
i = ~(J-Ll)xl-. 


PJ . A J . 
~l5L JXS = np S 3XSl 

~^LS 3Xl = ®‘ p lS 3Xl ' 


( 8 ) 

(9) 


Therefore the angular frequency with which s and l pre- 
cess around the total angular momentum j is given by 


n — ii p ^ — fi p 



cos<^+ 1 


= % 



- COS 4> 


+ sin 2 cj). 


J 2 M s R 2 s (s-l) 


2 


( 10 ) 
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For a circular orbit (e = 0), Cl reduces to 


D = —ncosd—Jo — 



sin 2 </>, (11) 


with 

L 1 M p n / a \ 2 
S = GM~ s CJ~ s \R s ) ’ 


( 12 ) 


where C is the moment of inertia coefficient of the central 
star and w s = 2tt/P s is the stellar spin frequency (Masuda, 
2015). It can be shown that Cl of equation (11) is identical 
to the precession frequency in B13; see their equations (6) 
or (8). 

The standard stellar evolution theory suggests that T- 
Tauri stars do not yet develop an internal radiative core. 
Thus PTFO 8-8695 is expected to be fully convective, and 
J 2 in this case is given by (Murray & Dermott, 1999) 


Ji 


Cf = C 


2GM S 


2 \ n ) \ a 


where 


/ 


-*^s,eq 


s,pol 


R* 


, ; 2 

1 b s,eq 

2 GM S 


(13) 


(14) 


is the oblateness of the central star, and i? s ,e q and i? SjPO i 
are stellar equatorial and polar radii, respectively. 

Since Kepler’s third law is written as 


Rs = f 3 n 2 \ 1/3 
a \4:TrGp s ) 


(15) 


f l in equation (11) is determined by the six parameters 
(n, (f>, oj s , p B , C , M p /M s ). Among these, the first four are 
constrained from the transit shape, and we adopt C = 
0.059, the value for the Sun, following B13. Thus M p /M s 
is the only parameter that is additionally constrained from 
the precession analysis. 

Note that since the equations of motion depend on p s 
alone, we cannot evaluate the stellar mass (M s ) and ra¬ 
dius (Rg) separately. We also emphasize here that the 
current analysis does not fix the absolute mass or radius 
( M s ,M p ,R s ,Rp ), but their ratios (M p /M s ,R p /R s ) only. 


4. Data analysis: methods and results 

4-1. Data reduction 

The photometric data analysed by B13 are kindly pro¬ 
vided to us by Julian van Eyken, which include 11 reliable 
transits in 2009 and 6 in 2010 observed at 0.658 pm. The 
effects of stellar activity, long-term periodical fluctuations 
(for instance, due to stellar spots) and background noise 
are already removed and whitened. 

We repeat the data reduction ourselves again following 
B13: 

1. Phase-fold. All transit lightcurves are stacked with 
the orbital period of 0.448413 days. 


[Vol. , 

2. Clipping. During orbital period of 0.448413 days, 
the transit actually continues at best for ~0.2 days. 
Thus a major part of the phase-folded lightcurve 
corresponds to out-of-transit phase. In order to re¬ 
duce the number of the data points and subsequent 
computational cost, we focus on ± 0.1 days around 
the transit center. 

3. Binning. All data points are combined into one- 
minute bins. 

4. Positioning the 2009 and 2010 lightcurves. B13 
found that the transit centers for 2009 and 2010 
data are 30861700 and 60848300 s measured from 
2009 January 1 UTC. Therefore we position 2009 
and 2010 lightcurves at those epochs. The mutual 
separation between the 2009 and 2010 lightcurves 
corresponds to 774 transits. 

4-2. Parameter fitting procedure 

Next we perform the parameter fitting as follows: 

1. Calculate Cl assuming a set of values for the param¬ 
eters listed in Table 1. 

2. Calculate s and l at the observed epochs in 2009 
and 2010 by extrapolating those at t c through the 
precession model prediction. 

3. Evaluate angular parameters (i s , * 0 rb, Cl) at the ob¬ 
served epochs in 2009 and 2010. 

4. Compute the model lightcurves at those epochs us¬ 
ing a routine developed by Masuda (2015), and 
evaluate y 2 = Xb(flux da ta,i -Auxmodei.i) 2 /^^ be¬ 
tween model and phase-folded lightcurves, where i 
indicates the variables evaluated at time t, in 2009 
and 2010. 

The above procedures are iterated so as to find the set 
of parameters that minimizes y 2 with the Levenberg- 
Marquardt algorithm (Markwardt, 2009). In what follows, 
we discuss in terms of the reduced chi-square y 2 (i.e., y 2 
normalized by the degree of freedom of the analysis of 
565). 

We see how the minimum y 2 value differs on P s — 
Mp/M s plane in our analysis. This strategy makes it easer 
to compare the results with and without the spin-orbit 
synchronous condition. The result of B13 corresponds to 
P s = 0.44841 days line. Also future constraints on M p from 
radial velocity measurements will be easily incorporated 
in this plane. 

We search for the best-fit parameters over the follow¬ 
ing regions: 0.12 < p s (g/cm 3 ) < 1.21, 0 < R p /R s < 0.5, 
7t/2 < i s < 7r, 0 < « or b < 7r, and 0 < fl < 2n. The lower and 
upper limits of p s are adopted from B13. Note that the 
symmetry in the gravity-darkened lightcurves introduces 
the degeneracy on the angular configuration of the system 
(Barnes et ah, 2011); (z s , ? or b, fi) <t=> (n — i B , tt — i or b, — Cl). 
While i s is originally defined between 0 and 7 r, therefore, 
we confine * s between n/2 and n. We abandon those so¬ 
lutions that violate the stellar break-up condition beyond 
which the star is not stable gravitationally: 
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P s > 2 tt 


4nGp s 


^.3. Results 

Figure 2 shows the minimum Xr values on the P s - 
M p /M s plane. The upper panel shows the three dimen¬ 
sional view, while the lower panel illustrates the corre¬ 
sponding plot as a color-coded contour map. The analysis 
of B13 was performed along the line of P s = 0.44841 days, 
and their best solutions are plotted as a yellow circle for 
M s = 0.34 M 0 and a yellow triangle for M s = 0.44 Mq. 

Figure 2 suggests that slightly better, or at least equally 
good, fits are obtained outside the synchronous condition 
(. P s = 0.448413 days), i.e., in the more rapid stellar spin 
and massive planet regime. Therefore, the synchronous 
condition over-constrains the system parameters, at least 
unless we have more precise observational data for the 
stellar spin period, P s . 




0.4 0.5 0.6 

P (days) 


P s or Mp/Mg. Now one can easily recognize that Xr is 
very sensitive to Pq, and the acceptable solutions are lo¬ 
cated around Pq = 199, 475, and 827 days (see the lower 
panel in Figure 3). For simplicity, we refer to them as 
short, middle, and long solutions, respectively, in the fol¬ 
lowing discussion. The best-fit parameters for the three 
solutions are summarized in Table 2 together with those 
of B13. Table 2 indicates that our solutions prefer the 
more rapid stellar spin and massive planet regime than 
those of B13. 



Pa (days) 



Fig. 3. Precession period dependence of \r- tied crosses are 
the data points on P s -M p /M s plane in Figure 2 arranged in 
terms of precession period. The lower panels are the zoom-ins 
around P^=199, 475, and 837 days in the upper panel. 


Fig. 2. Minimum Xr values for various values of P s and 
M p /M s . The upper panel shows the three dimensional distri¬ 
bution of minimum x'r- while lower panel is the correspond¬ 
ing color contour map. The results of Barnes et al. (2013) 
are those along the line of P s = 0.44841 days. Best solutions 
found by Barnes et al. (2013) are denoted as a yellow circle 
for M a = 0.34Mq and a yellow triangle for M a = 0.44 Mq. 

It is difficult to see from Figure 2, but the resulting 
solutions are roughly divided into three groups with dif¬ 
ferent precession periods (Pq = 2n/Cl). In order to clarify 
the point, we plot Xr against Pq in Figure 3, instead of 


The reason why we have three different solutions can 
be understood from Figure 4, which plots the predicted 
lightcurves for the three solutions; the short, middle, and 
long solutions are plotted as red, green, and blue lines, 
respectively. Because the transit period is very short, the 
lightcurves in the lower panel basically show the evolu¬ 
tion of the transit depth for three solutions. The five ver¬ 
tical lines in the lower panel are drawn at three-month 
intervals, each corresponding to the enlarged graph at the 
upper panels. By construction, all the three solutions re¬ 
produce the 2009 and 2010 data very well (upper-left and 
upper-right panels), and the significant difference of the 








































































Table 2. Best fit parameters for our three solutions and those in Barnes et al. (2013). 


parameter 

short solution 

middle solution 

long solution 

B13 (M s =0.34M e ) 

B13 (M s =0.44M e ) 

Ps(g/cm 3 ) 

0.32 + 0.01 

0.27 + 0.01 

0.30 + 0.01 

0.43 + 0.01 

0.57 + 0.02 

P s (days) 

0.390 + 0.008 

0.367 + 0.006 

0.331 + 0.008 

0.448410 

0.448413 

i-(°) 

123.9 + 2.2 

128.2 + 1.7 

131.3 + 1.5 

119.4 + 0.3 

120.3 + 1.3 

M p /M s 

0.0129 + 0.0014 

0.0178 + 0.0009 

0.0199 + 0.0008 

0.0084 + 0.0006 

0.0078 + 0.0007 

Rp/R s 

0.169 + 0.003 

0.208 + 0.008 

0.289 + 0.016 

0.159 + 0.007 

0.164 + 0.007 

*orb(°) 

62.9 + 0.8 

55.3 + 0.5 

50.9 + 0.6 

65.2 + 1.6 

69.3 + 1.3 

n(°) 

133.4 + 3.3 

131.2 + 2.4 

130.2 + 2.2 

136.1 + 5.2 

125.5 + 0.5 

<K°) 

75.3 + 2.5 

85.8 + 1.8 

92.3 + 1.6 

69 + 3 

73.1 + 0.6 

p n( days) 

198.7 + 15.6 

474.6 + 21.1 

826.9 + 53.3 

292.6 

581.2 


NOTE. - Our angular parameters (i s ,i 0 rb,ty are evaluated at t c = 2455543.9402 HJD while those of B13 are at 
f c =2455536.7680 HJD, whose interval corresponds to 7.1722 days or 16 transits. Spin-orbit angle <fi is time-invariant. 


lightcurves shows up only during the unobserved epochs. 

The examples plotted in the three upper-middle panels 
illustrate the importance of frequent monitoring of the 
system in order to specify the system parameters pre¬ 
cisely. Indeed some phases without any transit signals 
are predicted during their whole evolution because plan¬ 
etary orbital inclination significantly deviates from 90°. 
Thus measurements of the transit depth, even with a rela¬ 
tively low signal-to-noise ratio, greatly help distinguishing 
among the three solution groups. 

4 - 4 - Preliminary comparison with photometry at 
Koyama Astronomical Observatory 

In this subsection, we attempt a preliminary compar¬ 
ison with the transit photometry taken with Araki tele¬ 
scope at Koyama Astronomical Observatory (KAO) on 
2014 November 23, 27, December 2, 23, and 2015 January 
10. The photometry has been performed by a dual-band 
imager, ADLER (Araki telescope DuaL-band imagER, 
2048 x 2048 pixels and 12 arcmin 2 FOV for each CCD cam¬ 
era). The observation was performed with a 670 nm-split 
dichroic mirror, and the SDSS g'- and i'-band filters. The 
exposure time is set to 180 s, and 2 exposures are stacked. 
The observational data are reduced with standard proce¬ 
dures for photometry, i.e., dark-subtraction, flat-fielding, 
and aperture photometry (daophot task in IRAF). For rel¬ 
ative photometry, we adopt a star located at (RA, Dec) 
= (05:25:07.89, +01:31:16.1) as a reference star. We also 
compared the results with other reference stars, and con¬ 
firmed that the lightcurve did not change significantly, less 
than 1%. Since g'-band data do not have enough signal 
for the current purpose, we only focus on i'-band data in 
this paper. 

The bottom panel in Figure 5 shows the predicted evo¬ 
lution of transit depth for the three solutions over the 
observed epoch. The five vertical black lines indicate the 
above five days, whose enlarged views are shown in the 
upper panel. The predicted transit depths at those dates 
are less than 1% for the short solution and greater than 
approximately 10% for the middle and long solutions. The 
filled circles, on the other hand, indicate the data points 
after removing the long-term stellar activity. Since there 


are no clear transit signals beyond the photometric noise 
level of a few percents, the data seem to prefer the short 
solution. Although this conclusion is admittedly very pre¬ 
mature, it clearly indicates that such fairly crude photom¬ 
etry may be very useful in distinguishing among the three 
solutions. We plan to conduct the photometric observa¬ 
tions again during winter in 2015 and 2016. 

5. Long-term dynamical stability 

The results and discussion so far are all focused on the 
short-term dynamics of the system; the nodal precession 
whose expected time scale is several hundred days. In ad¬ 
dition to such a short-term behavior, long-term dynam¬ 
ics also contributes to getting an insight into the nature 
of the pre-main-sequence star and close-in planet system. 
When pursuing the longer-term evolution, tidal effect be¬ 
tween the star and close-in planet should be considered. 
In the PTFO 8-8695 system, particularly, the star and 
planet are expected to be strongly influenced by the tidal 
effects because the planetary semi-major axis is very small 
(a < 2R S ). 

We investigate the long-term stability of this system 
following the equilibrium tidal model by Correia et al. 
(2011). Note here that the two new parameters show up 
to describe the efficiency of the tidal effect; Love number 
and tidal delay time At, both of which refer to the 
fluid property of the central star. Tidal evolution of the 
system is computed by numerical integration of equations 
(19) - (21) in Correia et al. (2011). The efficiency of tidal 
effect is also represented by the reduced tidal quality fac¬ 
tor Q , which is written in terms of the above parameters 
as Q= 4raj . 3 2 At (Barker & Ogilvie 2009, Lai 2012, Rogers 
& Lin 2013, Xue et al. 2014). We adopt our short solu¬ 
tion as initial conditions and the values of Love number 
and tidal delay time estimated for Sun-like stars (&2 = 
0.028 and At = 0.1 s following Correia et al. 2011, which 
gives 10 6 for an orbital period of 0.448413 days). The 
upper, middle, and lower panels of Figure 6 show the re¬ 
sulting evolution of planetary semi-major axis, spin-orbit 
angle, and stellar spin and planetary orbital periods, re¬ 
spectively. Among them, it is worth noting that tidal ef- 
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Fig. 4. The evolution of the transit lightcurves for three 
possible solutions in Table 2 from 2009 to 2010 observa¬ 
tional epochs. Short, middle and long solutions correspond 
to red, green, and blue lines, respectively. Lower panel il¬ 
lustrates the whole evolution from 2009 (leftmost vertical 
line) to 2010 (rightmost vertical line). The five vertical lines 
are drawn for about every three months, each correspond¬ 
ing to the enlarged graph at the upper panel. Enlarged 
graphs for 2009 (leftmost) and 2010 (rightmost) are with 
the observational data (black points), while the interme¬ 
diate three are only theoretical calculation of lightcurves. 

feet leads to the planetary orbital decay (a —>• 0) with the 
time scale less than 10 4 yrs (upper panel). This time scale 
is smaller than the age of PTFO 8-8695 system (2.6-2.7 
Myrs) by at least two orders of magnitude. 

This rapid orbital decay is totally inconsistent with the 
observational picture that favors the survival of a close- 
in planet in the PTFO 8-8695 system. In addition, both 
spin-orbit angle </> and orbital period P 0 rb evolve quickly 
in accordance with the rapid orbital decay (middle panel 
and blue line in lower panel), which is also incompatible 
with the significant spin-orbit misalignment or spin-orbit 
near synchronization obtained by the lightcurve analysis. 
Since &2 = 0.028 and At = 0.1 s are estimated for Sun- 
likc stars, these issues suggest that the internal structure 
or fluid property of PTFO 8-8695 is significantly different 
from those of the main-sequence stars. Reasonable con- 


Fig. 5. Same as Figure 4, but for the observational 
term of KAO from 2014 November to 2015 January. 



Fig. 6. Tidal evolution of the planetary orbital semi-major 
axis (top), spin-orbit angle </> (middle), and stellar rotational 
and planetary orbital period (bottom). In the upper panel, 
stellar inside ( a/R s < 1.0) is described. In the lower panel, 
stellar rotational and planetary orbital periods are shown 
as red and blue lines, respectively. This calculation em¬ 
ploys the Love number of &2 = 0.028 and tidal delay time 
of At = 0.1 s, both of which are estimated for Sun-like stars. 
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straints on k 2 and At of PTFO 8-8695, therefore, allow us 
to approach the internal structure and fluid properties of 
pre-main-sequence stars, which are poorly known today 
due to the scarcity of the observed samples. 

6. Summary and discussion 

The spin-orbit angle <f> plays a key role in investigat¬ 
ing the formation mechanism of hot Jupiters (e.g., Crida 
& Batygin 2014 and references therein), and then var¬ 
ious methods to determine the angular configuration of 
the planetary system are developed and utilized. The 
most popular one is to use the Rossiter-McLaughlin (RM) 
effect (Rossiter 1924, McLaughlin 1924), which reveals 
the sky-projected spin-orbit angle through spectroscopy 
(e.g., Ohta et al. 2005, Hirano et al. 2012). By combining 
the RM effect and asteroseismology, furthermore, three- 
dimensional spin-orbit angle <f> can be determined (Lund 
et al. 2014, Benomar et al. 2014). 

The stellar gravity-darkening effect provides another 
methodology to measure <j> through the analysis of photo¬ 
metric transit lightuchrves for rapid rotators (Barnes et al. 
2011, Zhou & Huang 2013, Aiders et al. 2014, Masuda 
2015). Advantages of this method are that one can de¬ 
termine three-dimensional (i.e., not projected) spin-orbit 
angle directly and that the analysis is possible without the 
spectroscopical data (i.e., all we need is the photomet¬ 
ric transit lightcurves). When nodal precession occurs, 
furthermore, more precise constraints on the system are 
obtained (Philippov & Rafikov 2013, Barnes et al. 2013, 
Masuda 2015). 

Being a rapid rotator, PTFO 8-8695 is therefore an 
ideal sample to check the model of gravity darkening and 
nodal precession. We reanalyse the time-variable tran¬ 
sit lightcurves of this planetary system with the above 
theoretical model, discarding the spin-orbit synchronous 
condition assumed in the previous analysis, Barnes et al. 
(2013). We find three different groups of solutions that re¬ 
produce the observed data, whose precession periods are 
199 ± 16, 475 ±21, and 827 ±53 days. These solutions are 
slightly better, or at least equally good compared to those 
in the previous work. The three solutions reproduce 2009 
and 2010 observational photometry very well, whereas the 
theoretically-predicted transit lightcurves at unobserved 
epochs are found to be totally different. Difference in 
transit depth is particularly useful to distinguish among 
them with the frequent monitoring of the system. We 
also present a preliminary comparison of the lightcurve 
prediction with the additional photometry taken with 
Araki telescope at Koyama Astronomical Observatory, 
and find that the solution with the precession period of 
199 ± 16 days seems to be preferred. The solution im¬ 
plies p s = 0.32 ± 0.01 g/cm , P s = 0.390 ± 0.008 days, 
M p /M s = 0.0129 ± 0.0014, R p /R s = 0.169 ± 0.003, and 
(f> = 75.3°±2.5°. Although this result is still preliminary, 
more precise constraints on the system configuration are 
expected to result from our continuous photometric ob¬ 
servations during winter in 2015 and 2016. 

The long-term analysis taking into account the tidal 


[Vol. , 

effect between star and planet in the PTFO 8-8695 sys¬ 
tem suggests that tidal parameters &2 and At for pre- 
main-sequence stars are incompatible with those for main- 
sequence stars; otherwise, the orbiting planet is engulfed 
by the central star on the very short time scale. This fact 
implies that the internal structure of young stars is likely 
to be different from that for matured stars, and requires 
further investigation. 

Indeed, gravity darkening becomes more manifesting as 
the star rotates more rapidly. This is in contrast to the 
RM effect, which becomes harder to measure in the rapid 
rotators. Therefore gravity darkening provides a comple¬ 
mentary methodology to the RM effect in the measure¬ 
ment of </>. Since younger stars are known to rotate more 
rapidly than those matured, gravity darkening is available 
to unveil the distribution of (f> for the younger systems 
by the future observations. Although tidal evolution of 
<j> causes the initial memory on the formation process of 
hot Jupiters to be lost (Barker & Ogilvie 2009, Lai 2012, 
Xue et al. 2014), one can minimize this undesirable effect 
by preferentially selecting younger systems as targets to 
measure <j>. These aspects suggest that gravity-darkening 
(combined with nodal precession, if possible) model is very 
useful to study the formation scenario of close-in planets. 
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